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Emitting Lasers 



On 
On 



(N 
(N 



O 

• i-H ■ 

-f— > ■ 

Or 

q ; 

o ■ 
^ ■ 

Oh: 



> 

o 
in 
o 

5 ■ 

o : 

^ : 

Ph. 



X 



Holgcr F. Hofmann and Ortwin Hess 
Institute of Technical Physics, DLR 
Pfaffenwaldring 38-40, D-70569 Stuttgart, Germany 
(February 2, 2008) 



We present a generalized rate equation model for the array 
and polarization dynamics of general one or two dimensional 
arrays of Vertical Cavity Surface Emitting Lasers (VCSELs) . 
It is demonstrated that our model includes both the previous 
theory for edge emitting laser arrays and the theory of polar- 
ization dynamics in quantum well VCSELs in a single unified 
description. The model is based on the physical assumption 
of separated carrier densitiy pools individually coupled to dif- 
ferent light field modes. These modes, defined by the sepa- 
rate gain profile of each pool, interact through the coherent 
dynamics of the light field derived from Maxwell's equations 
and represented by the coefficients for index and loss guiding. 
The special case of two densities and two light field modes 
is solved and the implications of the results for large VCSEL 
arrays are discussed. Our analytic results show that typical 
solutions of the split density model range from phase locking 
to chaos, depending on the magnitude of the frequency coef- 
ficients. For weak coupling, the stable supermode is always 
the mode of highest frequency. This indicates that anti-phase 
locking is the only stable phase locking possible in arrays. 



I. INTRODUCTION 

If vertical cavity surface emitting lasers (VCSELs) are 
fabricated in the form of a densely packed two dimen- 
sional array as shown in Fig.^, their highly symmetric 
geometry around the axis of emission calls for an inves- 
tigation of both polarization dynamics and laser array 
interactions. Previously these two aspects have been 
treated separately. Experimental investigations of VC- 
SEL arrays |l]-|| usually cite theories derived for the one 
dimensional case of side-by-side edge emitting laser ar- 
rays such as the theory introduced by Winful and Wang 
in 1988 j|. On the other hand, a model for the descrip- 
tion of polarization dynamics in individual VCSELs was 
introduced by San Miguel, Feng and Moloney in 1995 ||. 

Both models are rate equation models for type B lasers. 
Adiabatic elimination of the dipole density is used to re- 
duce the description of the electronic system to carrier 
densities. However, the similarity does not end here. 
In both models, the carrier density is split into sepa- 
rate pools, each of which couples to a different light field 
mode. The differences between the models arise only 
from the choice of parameters which describe the dynam- 
ics of the electromagnetic field and the diffusive exchange 
of carriers between the carrier density pools. 



Since the physical assumptions of the models for polar- 
ization and array dynamics are so similar, we will show in 
the following that indeed a general model for the polariza- 
tion and array dynamics of arbitrary VCSEL arrays can 
be formulated. The coefficients of this generalized split 
density model will then contain all the information nec- 
essary to describe any special case. Therefore, the choice 
of the coefficient, which describe the physical properties 
of the device studied, is critical and a good understand- 
ing of the microscopic physics from which they can be 
derived is necessary to distinguish between common and 
exotic cases. 

In section ||, the physical justifications for introducing 
split densities in both the polarization and the array dy- 
namics is discussed. Parameters for carrier diffusion and 
scattering are approximated. 



In section II] 



we present 

the general formalism of the split density model. The 
coefficients derived from MaxwelPs_equations are intro- 
duced and discussed 



In section IV 



the special case of 
two densities and two light field modes is treated. Fixed 
points and limit cycles are analytically calculated. In sec- 
tion [v|, our analytic results of the two density model are 
applied to large arrays. A stability analysis is presented 
and conditions for chaotic dynamics are discussed. In 



section VI, we conclude by discussing the strengths and 



weaknesses of the split density model as the basis of a de- 
scription for general semiconductor laser array dynamics. 



II. THE PHYSICAL JUSTIFICATION OF SPLIT 
DENSITIES 

A. Polarization dynamics: separation of carrier 
density pools by angular momentum 

In quantum well VCSELs, the axis of symmetry of 
the quantum well structure is identical to the axis of 
laser light emission. Therefore, the conservation of an- 
gular momentum may be applied to describe the emis- 
sion process. Photons with an angular momentum of +1 
(right circular polarized light) may only be emitted by 
electron-hole pairs with a total angular momentum of 
+1 and photons with an angular momentum of -1 (left 
circular polarized light) may only be emitted by elec- 
tron hole pairs with a total angular momentum of -1. If 
only the bands closest to the fundamental gap are con- 
sidered, electrons have an angular momentum of ±1/2 
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and holes have an angular momentum of ±3/2. Conse- 
quently, electrons with an angular momentum of +1/2 
may only recombine with holes of -3/2 angular momen- 
tum, emitting left circular polarized light. Electrons with 
an angular momentum of -1/2 may only recombine with 
holes of +3/2 angular momentum, emitting right circular 
polarized light. 

Due to the conservation of angular momentum, the 
carrier density is therefore split into a carrier density 
emitting only right circular polarized light and a car- 
rier density emitting only left circular polarized light. 
This effect was first included in a rate equation model 
by San Miguel and coworkers in ||. In that model, it is 
assumed that the two densities interact by spin flip scat- 
tering processes which give rise to a diffusive exchange 
of carriers between the two densities. Experimental in- 
vestigations of quantum well structures indicate that this 
scattering rate is of the order of 10 10 s _1 , similar to the 
carrier recombination rate The assumption of split 
densities should therefore be well justified if the separa- 
tion of bands in the quantum well structure is sufficient 
to exclude the participation of higher bands in the laser 
process. 



approximation of the real physical process of continuous 
diffusion which is quite different from the spin flip scat- 
tering involved in the carrier exchange between carrier 
density pools of different angular momentum. However, 
as long as the light field is only sensitive to the average 
carrier density within one laser of the array and higher 
transversal modes are neglected, the overall effect of dif- 
fusion can be summarized in a single carrier exchange 
rate. 

We cannot but note that some of the coupled mode the- 
ories frequently cited use the assumption of homogeonous 
and stationary gain, ignoring the dynamical effects of 
spatially resolved carrier densities altogether Jj], 3|,p|-fio]. 
In type B lasers, that assumption is only valid at the 
fixed points. As a consequence, those theories can not 
be used as starting points for an analysis of the stability 
and the dynamics of laser arrays. 



III. GENERAL FORMALISM FOR THE 
DESCRIPTION OF VCSEL ARRAYS 

A. Rate equations for the carrier densities: 
injection, non-radiative recombination and diffusion 



B. Array dynamics: spatial separation of carrier 
density pools 

In a typical semiconductor laser array, carriers are 
injected into spatially separated regions of the active 
medium. If the distance of separation is sufficiently large, 
gain guiding effects will give rise to spatially separated 
optical modes, each associated by the gain function with 
its own carrier density pool. The localized modes are 
then coupled by Maxwell's equations, since the evanes- 
cent parts of the neighbouring optical modes overlap. 

The rate equation theory for laser arrays of this type 
was developed by Winful and Wang in 1988 [Q. However, 
in that early formulation, carrier diffusion is completely 
neglected. To check the justification of that assump- 
tion, we need to estimate the time it takes carriers to 
diffuse from one localized carrier density pool to a neigh- 
bouring one. Noting that the microscopically computed 
ambipolar diffusion constant of charge carriers in GaAs 
based lasers is typically about l/um 2 /ns = 10 cm 2 /s 0, 
we can estimate the diffusive carrier exchange rate be- 
tween lasers separated by a few micrometers to be close 
to 10 10 s _1 , similar to the spin relaxation rates expected 
for the polarization dynamics. However, the carrier ex- 
change rates depend on the square of the distance be- 
tween the lasers, which implies that the split density 
model rapidly loses its validity as the separation between 
elements of the array is reduced. In order to investigate 
the effects of diffusion, the possibility of carrier exchange 
between the densities is inc lud ed in the generalized split 



density model presented in III 



In this model, carrier diffusion is described in terms of 
a carrier exchange rate. Of course, this is just a crude 



The carrier densities of a VCSEL array can be sepa- 
rated into densities Di associated with a position in the 
array and either left or right circular polarization. In a 
two dimensional VCSEL array, i represents the coordi- 
nates (r, c, ±) for row number r, column number c and 
plus for right circular polarization, minus for left circular 
polarization. When no electromagnetic field is present, 
the carrier dynamics including the effects of injection, 
non-radiative recombination and diffusion are 



dt 



-Di 



(1) 



where fii is the effective local injection current, ji the 
local non-radiative recombination rate and 2y = Tjj is 
the rate of diffusive carrier exchange between the den- 
sities Di and Dj. Note that the model assumes trans- 
parency at fi = 0. When comparing the predictions of 
the model with experiment, it must be taken into ac- 
count that /i represents the additional injection current 
after transparency has been reached. The real injection 
current is /i plus a device dependent constant. The non- 
radiative recombination rate and the carrier diffusion can 
conveniently be combined into a single relaxation matrix 
of the densities by defining 



Tij = S,j (7, + ) T ik ) - Ti 



k 



(2) 



Equation |2J describes the diffusive relaxation to an 
equilibrium of injection and recombination of electron- 
hole pairs. This constitutes the complete carrier dynam- 
ics at zero field. Typical relaxation rates are around 
10 10 s _1 , as was discussed in [Il| 
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B. Maxwell's equations of the electromagnetic field: 
gain, loss and coherent propagation 

The electromagnetic field inside the optical cavity of 
the laser array can be described by the complex field 
amplitudes of localized modes Ei of either left or right 
circular polarization. Each mode is associated with one 
density pool. In a two dimensional array, i represents the 
same coordinates (r, c, ±) given above. This is an essen- 
tial feature of the split density model. The split densi- 
ties correspond to a set of modes for the electromagnetic 
field which is different from the Eigenmodes of Maxwell's 
equations. Instead, these modes are Eigenmodes of the 
gain. They are the natural modes for the description of 
relaxation oscillations and limit cycles as well as the sta- 
bility of fixed points. We will therefore chose this set 
of Eigenmodes and not the supermodes to describe the 
laser array. We note that we limit our theory to the fun- 
damental mode of each laser in the array and disregard 
possible effects of spatial hole-burning within each laser. 
A generalization to include higher modes is possible, but 
would make it more difficult to distinguish array effects 
from the internal dynamics of individual lasers. 

The dynamics of the electromagnetic field, includ- 
ing gain, loss and coherent propagation are given by 
Maxwell's equations, which are linear wave equations. In 
a general mode representation, they may be written as 

—Ei = WiD t (l - ia)Ei - + iu>ij)Ej. (3) 

j 

The carrier density dependent gain is given by WiDi. 
The rate of spontaneous emission into the laser mode Wi 
is typically around 10 6 s _1 . Since most arrays will have 
very similar values for the different Wi, the average and 
the deviations from this average may be seperated into 
Wi = w(l + fi). In this context, /, is the measure of the 
relative fluctuation in the coupling between density pool 
i and mode i. 

The alpha factor a describes band structure effects 
which shift the gain profile to higher frequencies as the 
total density is increased. Typical values range from 2 to 
6. Effectively, this shift of the gain profile introduces a 
carrier density dependent index of refraction. 

Kij describes the losses through the mirrors. In most 
cases, the coefficients coupling different modes are very 
small and most of the losses are given by K&y with R 
around 10 s . Small irregularities and birefringence in 
the mirrors may cause non negligible loss guiding effects 
which can be represented by small corrections Rlij. 

The matrix u>ij describes the unitary coherent dynam- 
ics of Maxwell's equations based on the index of refrac- 
tion of the array structure. Neighbouring modes interact 
by their overlapping evanescent waves, while right and 
left circular polarization modes are coupled by stress in- 
duced birefringence commonly found in epitaxially grown 
semiconductor heterostructures |l5| . Since the magni- 
tude of these effects is highly dependent on device archi- 



tecture, it is difficult to estimate. However, their order 
of magnitude should be roughly between 10 10 s _1 and 
10 12 s _1 to have a significant effect on the dynamics. For 
a given architecture, numerical simulations may be used 
to determine the effective coupling constant. 

For stationary densities Di , the dynamics of the elec- 
tromagnetic field is strictly linear. It is therefore possible 
to find orthogonal Eigenmodes of the dynamics includ- 
ing the cooperative effects of gain guiding (wiDi), loss 
guiding (k^) and index guiding (u>ij). These modes are 
the supermodes of the coupled mode theories. Note that 
Eigenmodes of different density distributions Di need not 
be orthogonal. Since semiconductor lasers are type B 
lasers, the dynamics of the Di cannot be neglected and 
the modes defined by stationary Di are next to useless 
for the investigation of dynamical effects and stability 
analysis. 

C. Combined dynamics of the carrier density and 
the electromagnetic field 

Now, the total dynamics of carrier density and elec- 
tromagnetic field can be constructed by adding only the 
term of induced emission, i.e. —2wiE* EiDi, to the car- 
rier density dynamics. The resulting equations of motion 
for N split densities and N complex field amplitudes rep- 
resenting, for example, an array of N/2 VCSELs with two 
polarizations each are 

j t Di = T v D l - + fi) E * E i D i (4a) 

j 

4$ = Ml + fi)Di(l - ia) - R)Ei 
at 

- y^X^kj + iu>ij)Ej. (4b) 

3 

These two lines of equations already contain a large num- 
ber of different physical effects. The localized modes in 
the array are not only coupled by overlapping evanescent 
waves (coij) but also by carrier diffusion (I\j) and by 
variations in the losses (hj). The injection current may 
be varied (/ii) and the coupling of the active medium to 
the laser mode may also be different for each laser in the 
array (fi). 

In the following, we shall focus our attention on ar- 
rays pumped homogeneously (/ij — fj.) and with negligible 
variations in the gain and loss coefficients (fi = = 0), 
resulting in 

jDi =»-J2 V ij D 3 - ^E*EiDi (5a) 

3 

^Ei = (wDi(l - ia) - R)Ei - iujijEj. (5b) 

3 

This equation describes the counteractive effects of co- 
herent field dynamics represented by u>ij and carrier dif- 
fusion represented by . 
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For the case of VCSEL polarization, the joint ef- 
fects of coherent coupling and non-negligible gain-loss 
anisotropy, kj ^ 0, has been investigated by the authors 
in a separate paper |1 1|. 



IV. THE TWO DENSITY MODEL 

A. Formulation of the rate equations 

To understand how the effects of coherent field dynam- 
ics and carrier diffusion influence the stability of phase 
locking and the transition to chaos, we will now start our 
discussion with an investigation of the most simple case 
of two carrier density pools and two light field modes. 
This case describes either the polarization dynamics of 
a single VCSEL or an array of two edge emitting lasers. 
In these contexts, the equations of the two density model 
have been studied previously by Martin- Regalado (polar- 
ization) |1| and by Winful and Wang (array) Q on the 
basis of numerical simulations. In the following we will 
review the results of the two density model and analyti- 
cally derive general properties of array and polarization 
dynamics. 

The equations for a symmetrical system of two carrier 
densitiy pools are defined by the parameters w, k, Th, 
Ti2, and uji 2 . It is useful to keep in mind that Tn is the 
sum of the carrier diffusion rate Ti2 and the non-radiative 
recombination rate 7. They read 



d 


= /x-r u D! + 


T12D2 - 


2wE^E 1 D 1 


(6a) 


d 

dt D > 


= /x-r u r> 2 + 


T 12 D l - 


2wE^E 2 D 2 


(6b) 


d F 
Jt El 


= {wDi{\ — ia 


) - k)Ex 


- lLUl 2 E 2 


(6c) 


d F 


= (wD 2 (l — ici 


) - fi)E 2 


- iu 12 E x . 


(6d) 



All rates are given by positive and real values. This is 
not a natural assumption for oj\ 2i since a phase trans- 
formation in one of the two modes would generally also 
change the phase factor of 0^12. A real and positive lo\ 2 
implies that the E\ +E 2 mode has a lower frequency than 
the Ei — E 2 mode for Di = D 2 (no gain effects on the 
mode structure). In an array, this is a natural choice. 
If the modes describe orthogonal polarizations, however, 
the phase factor of uji 2 is an important definition of the 
phase relation between the modes. 



B. Stability analysis 

If the two subsystems were not coupled by u>i 2 , they 
could be described as two identical lasers operating at 
the same injection current. Therefore, D\ = D 2 and 
El Ei = E2E2 is stationary. The introduction of lui 2 



removes the arbitrary choice of phase and defines the 
supermodes with E\ = —E 2 or Ex = E 2 as the only 
possible stationary solutions for the light field. 

The exact stationary solutions are given by D\ = D 2 = 
k/w and E\E\ — E^E 2 — /x/2k — ^/2w. To analyse the 
stability of the two possible supermodes, small fluctua- 
tions around these stable values are defined: 



SD = Di + D 2 - — 

w 

8n = El Ex + E*E 2 - 

d = D l -D 2 
S\ = El E\ — E^E 2 
s 2 = i{E* 1 E 2 -E* 2 E 1 ). 



R 



2 

w 



(7a) 

(7b) 

(7c) 
(7d) 
(7e) 



Note that in the case of polarization dynamics, si and 
s 2 are two of the three Stokes parameters. The third 
Stokes parameter is defined as S3 = E\E 2 + E^E\ and 
is equal to plus or minus the total intensity for the two 
stationary solutions. 

Using these parameters, the linearized dynamics 
around the stationary point of the antisymmetric super- 
mode with Ei — ~E 2 and S3 = — n s is 

( 5D \ 

Sn 



d_ 



(I 
Si 

S2 



\s 2 J 



1 -^ — wn s 


-2k 











wn s 

















-7- 


-wn s 2ri2 


-2R 











wn a 





W12 


V 





a.wn a 


-W12 












on 




d 




Sl 


) 


w 



(8) 

where n s = E^E\ + E^E 2 is the total intensity of the 
stationary solution. To obtain the stability matrix for 
the symmetric supermode with E\ = E 2 and S3 = +n s , 
one only has to change the sign of ui 2 . 

The stability matrix can be separated into two blocks. 
One two by two block describes the relaxation oscillations 
of the total intensity. The remaining three by three block 
contains all the relevant information on the interaction 
between the subsystems. For weak coupling, si and d 
show relaxation oscillations damped by diffusion, while 
S2 will be stabilized or destabilized exponentially by its 
interaction with d via a. 

The condition for the stability of a solution is, that the 
real parts of all Eigenvectors of the matrix given above 
are negative. For the symmetric supermode Ei + E 2 , 
which is the lowest frequency mode, this condition may 
be found by searching for a parameter set with one Eigen- 
value equal to 0. This is done by setting the determinant 
of the stability matrix to 0. The result of this analysis is 



lui 2 > 2na 



7 + wn s + 2T 



12 



(9) 
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Since we must assume that the diffusion rate Ti2 is not 
much larger than the non-radiative recombination rate 7 
for the split density model to be valid and since the rate 
of induced emission wn a cannot be much larger than 7 
for reasonable injection currents, this condition can only 
be fulfilled if 0J12 is as large as the loss rate R. This 
would require frequency separations of THz between the 
supermodes, which is unlikely if any reasonable gain, in- 
dex or loss guiding exists to separate the modes. Since 
such guiding effects are absolutely necessary to define 
the array, the low frequency solution is only stable if the 
description of the device as an array of coupled lasers 
breaks down. It is therefore reasonable to assume that in 
a laser array guiding effects reduce W12 to values below 
R. Consequently only the highest frequency supermode 
can be stable in a two mode array. If a different super- 
mode is stabilized, this is most likely not an array effect, 
but instead represents higher mode effects typical of a 
broad area laser device (l^Jl^]. In the case of polariza- 
tion dynamics in VCSELs, a birefringence effect of THz 
is unrealistically high. Experiments indicate that u>i2 is 
more likely to be of the same order of magnitude as 7 
]l5fl . Therefore, the stable polarization corresponds al- 
ways to the one with the highest frequency. Of course, 
the stability condition for the low frequency mode can 
still be fulfilled very close to threshold. This gives rise 
to a bistability which was reported in Jl2| . However, 
according to our estimates, this bistability should only 
occur extremely close to threshhold, where noise effects 
dominate the laser dynamics. This bistability is therefore 
not likely to be of any practical relevance. 

The condition for stablility of the anti-symmetric su- 
permode Ei — E 2 is found by searching for parameters 
permitting the existence of purely imaginary Eigenvec- 
tors. This indicates the point at which the relaxation 
oscillations of s\ and d become undamped. The condi- 
tion for damped relaxation oscillations is found to be 



au\2 < 7 



wn s 



2Ti 



(10) 



This condition is always fulfilled if the damping effect 
of the diffusion rate 2Ti2 exceeds the undamping effect 
of the coherent coupling auji2- Effectively, a moderate 
amount of diffusion stabilizes phase locking and prevents 
the transition to limit cycles and chaos. If awi2 exceeds 
7 + 2ri2, a minimal intensity is necessary to stabilize 
phase locking. Possibly, this intensity is unrealistically 
high, so that injection currents far above threshold would 
be necessary to stabilize the dynamics. Otherwise, a 
transition from stability to a limit cycle will be observ- 
able as the injection current is lowered. 



C. Dynamics close to the fixed point: damping and 
undamping of relaxation oscillations 



In order to provide a foundation for the discussion of limit 
cycles and for the extrapolation of the results to larger 
arrays, we will now take a closer look at the dynamics 
of fluctuations around the fixed point. This is done by 
finding the Eigenvectors and Eigenvalues of the stabil- 
ity matrix. Although, in principle, approximations are 
not necessary to diagonalize the two by two and three by 
three blocks of the stability matrix it helps the interpre- 
tation of the physics if we make use of the fact that for 
semiconductor lasers, the cavity loss rate R is typically 
two orders of magnitude larger than the non-radiative 
recombination rate 7. As a result, the fluctuations of 
semiconductor lasers are dominated by fast relaxation 
oscillations of the frequency v — (2Rwn s 

) 1 / 2 . The sta- 
bility matrix can then be diagonalized by treating the 
damping term 7 + wn s + 2Ti2 and the coherent coupling 
W12 as small perturbations of the relaxation oscillations. 
In lowest order, the Eigenvalues Aj and the left and right 
Eigenvectors and are 

^1/2 = -5(7 + wn s ) ± iv 



a V2 = 73 { T 



>l/2 
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(11a) 



V J 



A, = 



-awi2 



a 3 = ( -a 1 



b, = 



1 4 




Vi7 



(lib) 



A4/5 = -5(7 + wn s + 2Ti2 - au)i 2 ) ± - 




a 4/5 = 7 =( T iJ^ 1 



b4/5 " T2 



(11c) 



The analytic results of the stability analysis show that 
both semiconductor laser arrays and birefringent VC- 
SELs tend to emit light in the mode of highest frequency 



In (11a) and (11c), solutions 1 and 4 refer to the up- 
per sign and solutions 2 and 5 refer to the lower sign, 
respectively. 
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The solution 1/2 describes the linearized relaxation os- 
cillations of the total intensity and carrier density, solu- 
tion 3 describes the stabilization of anti-phase locking 
and solution 4/5 describes the relaxation oscillations of 
the difference in the intensities and carrier densities. 

The solutions 3 and 4/5 describe the supermode stabil- 
ity we are interested in. We note that an important fea- 
ture of the Eigenvectors is their non-orthogonality. This 
is caused by the a factor. The a factor converts density 
fluctuations into phase fluctuations. An example to illus- 
trate this effect can be given by the relaxation dynamics 
of a fluctuation in the intensity difference si, described 
at time t — by the vector (0,0,0,1,0). As shown in 
Fig. ||, this intensity fluctuation induces a phase fluctua- 
tion proportional to a as it relaxes. 

The relaxation rates of the solution are influenced by 
this correlation between phase fluctuations and intensity 
fluctuations. In solution 3, the a factor converts the time 
derivative of the intensity difference, d/dt si — W12S2, into 
a time derivative of the phase, A3S2 = a d/dt si, causing 
stability or instability, depending on the mode. In solu- 
tion 4/5, the time derivative of the phase is likewise con- 
verted into a time derivative of the intensity difference, 
increasing or decreasing stability. If the phase difference 
is stabilized, the relaxation oscillations are destabilized 
and vice versa. Mathematically, this is necessary because 
the sum of all Eigenvalues must be equal to the trace 
of the stability matrix. 

Since the a factor increases frequency with increasing 
carrier densities, the equations show that only the super- 
mode with highest frequency is stabilized in the presence 
of relaxation oscillations. However, the stability analysis 
has shown that strong coupling may destabilize this so- 
lution as well. In the Eigenvalues of solution 4/5, this is 
expressed by the negative contribution of au>\2- 

To understand the undamping effect given by au>\2, let 
us recall that S\ and S2 oscillate in phase, with S2 = as\. 
Consequently, the matrix element u>\2 causes an undamp- 
ing of this oscillation by adding a term equivalent to 
+au>\2Si to the time derivative of s\. Effectively, we 
thus encounter the surprising property that the coher- 
ent dynamics of the light field pumps intensity from the 
mode at low intensity to the mode at high intensity. This 
destabilizing effect is countered by the carrier density re- 
laxation effects when the carrier density difference d is 
maximal and the light field intensity difference s\ is zero. 
In this way, non-radiative (7) and radiative (wn s ) carrier 
recombination and diffusion (2Fi2) stabilize the relax- 
ation oscillations by damping d, while coherent coupling 
(awi2) destabilize them by undamping s%. 



D. Low amplitude limit cycle 

If the relaxation oscillations are undamped, the result 
is a low amplitude limit cycle. For small amplitudes, this 
limit cycle still resembles relaxation oscillations with 



/ 2ji 

d(t) = A$\ — — cos(vt) 
V wn s 

Si(t) — A sin(ut) . 



(12a) 
(12b) 



The damping term is still linear, since the diffusion and 
the carrier recombination terms are linear over a wide 
range of carrier densities. The undamping term is more 
complicated. It results from the phase dynamics 



tan(cj)) 



. E\E2 — E2 Ei S2 
''ElE 2 ~+E*Ei = H 



(13) 



of the phase difference tfi between E± and i?2- For 
small fluctuations, S3 —n s and tan{<j)) <f>, so that 
S2 ~ —n s (j). If non-linear terms are included, the phase 
dynamics of the limit cycle is given by 



a 



Aasin{vt) . 



(14) 



Since the total intensity n s , the 1st Stokes parameter 
si and the ratio between the 2nd and 3rd Stokes param- 
eters (i.e. tan(4>)) are known, we can determine S2 using 



si 



si 



s 2 (t) = ^Jnl- Sl {t) 2 sin(4>(t)). 



(15) 



The 3rd order correction to the linear dynamics of s 2 
is then given by 



s 2 (t) = a Sl (t)(l - Sl (t) 2 



a 



6n 2 



+ ...). 



(16) 



For small s\/n s , the term of 3rd order in s\ reduces 
and flattens the maxima of S2- To determine the time 
averaged effect on the amplitude Aq, the contributions 
to the time derivative of Si in phase with the oscillations 
must be found. This is the Fourier component of sin(vt). 
Since the Fourier component of sin[yt) in sirt'{vt) is 3/4, 
the total time averaged damping is 

!/ - ^ ^ 3 A%3 + a 2 , . 

-(7 + wn s + 21 12 - auj 12 ) + -awi2^r — -p, — ■ (17) 
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Setting this term to 0, the amplitude of the limit cycle 
is found to be 



A = 2n t 



j_ (1 _ 7 + ffl». + ZT U )j (18) 



3 + a 



OtUJ\2 



This equation describes the transition from stability to 
a limit cycle possible for strong coupling and low diffu- 
sion. The intensity dependent contribution to the stabil- 
ity is the rate of induced emission, wn s . 

The amplitude of the limit cycle which has been calcu- 
lated numerically in the paper on laser arrays by Winful 
and Wang Q can now be determined analytically using 
this formula. That limit cycle was calculated at an in- 
jection current of 1.1 times threshold, while stability is 
reached at 2 times threshold for the parameters used. 
Fig|| shows the dependence of Aq on the injection rate 
for this choice of parameters. 



G 



V. GENERALIZATION OF THE TWO DENSITY 
MODEL RESULTS TO LARGER ARRAYS 

A. Effective next neighbour interaction 

Our analysis shows that in the case of stability, two 
interacting lasers will emit light 180 degree out of phase. 
This result can also be applied to larger arrays. Usually, 
the coupling constants only couple next neighbours. 
If all u>ij are real and positive, the stable mode will be the 
one in which next neighbours emit light 180 degree out 
of phase. Note that in the case of a VCSEL array, the 
two circular polarizations must be considered as separate 
members of the array, coupled by birefringence. In this 
sense, VCSEL polarization adds a third dimension to the 
otherwise two dimensional array. 

In most conventional geometries, it is very easy to find 
a single stable supermode in which all next neighbours 
emit exactly out of phase. For example, a square array 
of VCSELs will emit light of a single linear polarization 
in the highest frequency supermode. In the far field, this 
corresponds to four intensity maxima emitted into direc- 
tions tilted towards the corners of the array. This result 
has been observed experimentally JjJ , providing evidence 
that the split density model is a valid description for re- 
alistic devices. 

In unconventional geometries such as triangular arrays 
or spatially varying birefringence, there may be more 
than one stable supermode. The stationary solutions 
may be obtained by letting the field amplitudes Ei be 
equal to +Eq or —Eq. Among these modes, promising 
candidates for stability are always the arrangements with 
the highest possible number of phase changes between 
next neighbours. 

Once the stable modes have been identified using the 
rule that next neighbours try to achieve anti-phase lock- 
ing, the stability of this mode with regard to phase fluc- 
tuations and relaxation oscillations may be investigated. 
Since we can assume that the fastest time scale is that 
of the relaxation oscillations, the two effects may be sep- 
arated. The justification for this separation is the same 
as the one underlying the stability analysis of IV: The 
Eigenvectors and Eigenvalues of the stability matrix may 
be determined by treating u>ij and the effects of damping 
and diffusion as small perturbations of the fast relaxation 
oscillations. The 3N dimensional system of equations for 
the carrier densities and the real and imaginary parts of 
the fields of N lasers can then be separated into a system 
with 2N dimensions describing relaxation oscillations and 
one of N dimensions describing phase coupling. 

In the following, we limit the coupling terms to next 
neighbour interactions, allowing a complete stability 
analysis of large symmetric laser arrays. 



B. Stability of the anti-phase locked supermode 

First, we will investigate the N dimensional stability 
matrix which describes the phase coupling of the ar- 
ray. In the two density model, the phase fluctuation 
was described by S2 ~ —n s 4>, where the phase differ- 
ence <p — &4>i — S(f>2 represents a measure of the difference 
between the phase fluctuations in system 1 and system 
2. If 5<f>i denotes the phase fluctuation in system i, the 
two density stability matrix for phase fluctuations is 



d_ 

dt 



2 

au)i2 
2 



aui2 
2 

2 



(19) 



This result may be generalized to larger arrays by using 
( |l9| ) as next neighbour interaction of phase fluctuations 
in arbitrary arrays. The general equation for the dynam- 
ics of phase fluctuations is then given by 



-c 
dt 



E 



au>. 



(54>j - 6<f>i). 



(20) 



This equation describes the linear relaxation of phase 
fluctuations in an array of arbitrary size and geometry. 
The stability matrix is defined by the awy. Note, that 
the relaxation has the properties of a diffusion such as 
the one described by for the carrier densities as in- 
troduced in III. This indicates that local phase changes 



diffuse to the neighbouring systems until the whole array 
is again phase locked. For an array with lattice constant 
ao, the phase diffusion constant is approximately given by 
a^aujij. Thus, for micrometer scale arrays, typical phase 
diffusion constants are probably around lfj,m 2 /ns. This 
is the same order of magnitude as the carrier diffusion 
constant, even though the physical properties defining 
each are quite independent of each other. 



C. Anti-diffusion and the route to chaos 

The remaining 2N variables of the stability matrix de- 
scribe relaxation oscillations of the carrier densities and 
the light field intensities of the array. 

In the two density model, the relaxation oscillation am- 
plitude of both the total intensity and the intensity differ- 
ence are damped by the non-radiative and the radiative 
recombination of carriers, i.e. 7 + wn s . In addition, the 
oscillations of the intensity difference are damped by dif- 
fusion, 2Fi2 and undamped by the coherent interaction 
between the two systems, auo\2. Again the coherent in- 
teraction acts as a diffusion process. However, this time 
the diffusion constant is negative: the coherent coupling 
gives rise to anti-diffusion, increasing the intensity differ- 
ence between next neighbours. 

If we denote the local amplitude of the relaxation os- 
cillations in the subsystem i by the complex value Ai , the 
two density result may be generalized to arbitrary arrays 
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on the basis of the same principles applied in the previ- 
ous subsection to the phase fluctuations. The dynamics 
of the Ai averaged over several relaxation oscillations is 
then given by 



dt 



4((7i + wn s )Ai 



-){Ai 



(21) 



Note, that Py has been decomposed into the contribu- 
tions from non-radiative carrier recombination 7, and 
from diffusion TJy, because they contribute to different 
types of effects, 7$ + wn s causes an overall relaxation 
of the oscillations, while the diffusion and anti-diffusion 
terms Ty — aoJij/2 are sensitive only to amplitude dif- 
ferences between nearest neighbours. For 7i = 70 = 
constant, the relaxation effect stabilizes all forms of re- 
laxation oscillations equally well. The diffusion stabi- 
lizes especially the relaxation oscillations with an am- 
plitude which varies strongly between neighbouring sys- 
tems. Consequently, these relaxation oscillations are also 
undamped most effectively when Xy < awy/2. In a 
regular lattice, the least stable relaxation oscillation in 
the presence of such anti-diffusion is the one in which 
Ai = —Aj for next neighbours. In a square lattice with 
constant coherent interaction w nn and constant diffusion 
T nn between next neighbours, the stability condition for 
this oscillation is 



< 2Tij + 7(7^ + wn s ) 



(22) 



When aujij exceeds this limit, the relaxation oscilla- 
tions become undamped and the near field intensity os- 
cillates in a chessboard pattern. Since this implies that 
the intensities of all VCSELs emitting in the same phase 
oscillate in parallel, this oscillation should show up as a 
small intensity peak in the otherwise dark center of the 
far field pattern. 

For stronger anti-diffusion, more and more relaxation 
oscillations become undamped and the amplitude of the 
oscillations increases. The chessboard pattern of the os- 
cillations will then be washed out and non-linear effects 
shift the frequencies and add higher harmonics. For very 
strong coupling, the neighbouring VCSELs will switch on 
and off in an unpredictable chaotic fashion. 



VI. PRACTICAL IMPLICATIONS AND OPEN 
QUESTIONS 

The split density model is the most simple way of de- 
scribing the array and polarization dynamics of VCSELs. 
It is also applicable to the array dynamics of conventional 
edge emitting arrays. 

For the most interesting case of large symmetric ar- 
rays, the results of the split density model are strikingly 
simple and clear: if the separation between the lasers in 



the array does not break down because of rapid diffu- 
sion, next neighbours lock into anti-phase laser emission. 
Consequently, only the highest frequency supermode is 
stabilized. 

For VCSEL arrays, this behavior has indeed been ob- 
served experimentally in In Q, this effect is ex- 
plained as the result of higher losses between the VCSELs 
in the array. In the split density model, this corresponds 
to a coupling via the £y terms introduced in III. Such 



terms do indeed stabilize the highest order supermode to 
a certain degree. However, since the ?y are proportional 
to the weak intensities, while the wy terms are propor- 
tional to the fields between the VCSELs , this effect is 
likely to be much smaller than the coherent interaction. 
We therefore conclude that the explanation given in jlj 
for the experimentally observed anti-phase locking may 
be a misinterpretation caused by the neglect of carrier 
dynamics in many coupled mode theories p)-|lO|]. Those 
theories certainly ignore the stabilizing effects of carrier 
dynamics and coherent interaction, which cannot be ne- 
glected in type B lasers. 

Another aspect of the split density model, pioneered 
in the numerical work by Winful and Wang jj] , is the ex- 
istence of limit cycles and chaos in semiconductor laser 
arrays. We have calculated the stability limit includ- 
ing carrier diffusion effects and have shown that stability 
is not as unlikely as predicted in the original paper by 
Winful and Wang. Indeed, carrier diffusion may stabi- 
lize phase locking. 

A careful analysis of the fluctuations in the center of 
a far field pattern such as the one observed in could 
reveal whether the array is well stabilized by diffusion 
(weak fluctuations in the center), shows anti-diffusion 
(fluctuations have a maximum in the center) or is al- 
ready oscillating in a weak limit cycle (coherent light is 
emitted into the center of the far field). Possibly, the 
correct range of parameters necessary for limit cycles and 
weak chaos could be realized experimentally by varying 
the distance between the VCSELs to change diffusion ef- 
fects and by using different types of optical guiding. For 
larger separations, pure gain guiding may be sufficient, 
while loss guiding is better for very small arrays. 

In the same spirit, VCSEL polarization may be investi- 
gated, using the quantum well width to change scattering 
rates and stress to induce birefringence. Here, too, the 
frequency spectrum of the polarization orthogonal to the 
laser mode reveals a lot about the mechanism of stabi- 
lization in the VCSEL and may actually include coherent 
emissions from a limit cycle. Possibly, some of the VC- 
SELs classified as having only poor polarization stability 
in reality show a well stabilized high frequency limit cy- 
cle and it is the time average of this limit cycle which is 
observed in polarization measurements. 

Of course, when interpreting experimental results, the 
underlying assumptions made in the split density model 
must be kept in mind. If the diffusion or the spin flip 
scattering is too strong, the split densities must again be 
considered as a single carrier density and the supermodes 
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become higher modes of a single laser with many modes. 
Consequently, the split density model can be used to de- 
scribe the transition from split to unified density. This 
has been attempted for polarization dynamics in quan- 
tum well VCSELs by Martin- Regalado in |L2|. However, 
that approach is only valid for a narrow range of parame- 
ters and it may be debatable if the limitation to the low- 
est quantum well bands implied in J5|,[l2| is still justified 
in quantum wells with such a strong spin flip scattering. 

If the injection current is increased, higher order modes 
of the individual lasers in the array may also participate 
in the dynamics. However, to simplify the presentation 
we have disregarded multi-transverse modes of an indi- 
vidual laser. When considering that case, one has to 
assume that the dynamics is not governed by array ef- 
fects, but instead by phenomena typical of broad area 
laser dynamics. The transition between the two regimes 
can be seen in numerical simulations of gain guided edge 
emitting arrays Jl3|,[l4| . In order to simulate the effect in 
the framework of the split density model, higher modes 
must be assigned to each density. However, this rapidly 
reduces the simplicity which makes the model so attrac- 
tive. 

In conclusion, we have shown that the split density 
model provides a simple description of semiconductor 
laser arrays and is equally useful for the description of 
VCSEL polarization. The results clearly show that the 
mode with the highest frequency is stabilized. This is in 
agreement with the experimental results reported in 
and indicates that the handwaving explanation of the 
anti-phase locking given in |l| may indeed be wrong. We 
have seen that strong coupling may undamp the relax- 
ation oscillations, leading through limit cycles to chaos. 
However, carrier difusion and both radiative and non- 
radiative carrier recombinations stabilize the anti-phase 
locked mode, which makes the stable operation of a VC- 
SEL array more likely than the chaotic case. 
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FIG. 1. Schematic geometry of a VCSEL array. 

FIG. 2. Relaxation dynamics of the phase variable S2 and 
the intensity difference s\ induced by a fluctuation in the 
intensity difference s\ at time t — 0. The choice of pa- 
rameters is fully characterized by a = 3, awi2 = u/15 and 
(7 + wn s + 2T\2 — aw\%)l1 = u/5. This example illustrates 
clearly, that fluctuations in the intensity difference si induce 
strong fluctuations in the phase difference described by S2- 

FIG. 3. Dependence of the limit cycle amplitude Aq on the 
carrier injection rate. The choice of parameters corresponds 
to that of [4] : a = 5, Ti2 = and 00)12 = 27. The injec- 
tion rate in units of the threshold injection rate is equal to 
1 + wn„/y. For this choice of parameters, the transition from 
the limit cycle to stabilization by induced emission is within 
the experimentally accessible range of injection current and 
would therefore be observable. 
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